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ABSTRACT 

Approximately 10-20% of all Kuiper belt objects (KBOs) occupy mean-motion resonances with 
Neptune. This dynamical configuration likely resulted from resonance capture as Neptune migrated 
outward during the late stages of planet formation. The details of Neptune's planetesimal-driven 
migration, including its radial extent and the concurrent eccentricity evolution of the planet, are the 
subject of considerable debate. Two qualitatively different proposals for resonance capture have been 
proposed — migration-induced capture driven by smooth outward evolution of Neptune's orbit and 
chaotic capture driven by damping of the planet's eccentricity near its current semi-major axis. We 
demonstrate that the distribution of comparable-mass, wide-separation binaries occupying resonant 
orbits can differentiate between these two scenarios. If migration-induced capture occurred, this 
fraction records information about the formation locations of different populations of KBOs. Chaotic 
capture, in contrast, randomizes the orbits of bodies as they are placed in resonance. In particular, 
if KBO binaries formed by dynamical capture in a protoplanetary disk with a surface mass density 
typical of observed extrasolar disks, then migration-induced capture produces the following signatures. 
The 2:1 resonance should contain a dynamically cold component, with inclinations less than 5-10°, 
having a binary fraction comparable to that among cold classical KBOs. If the 3:2 resonance also hosts 
a cold component, its binary fraction should be 20-30% lower than in the cold classical belt. Among 
cold 2:1 (and if present 3:2) KBOs, objects with eccentricities e < 0.2 should have a binary fraction 
~20% larger than those with e > 0.2. Other binary formation scenarios and disk surface density 
profiles can generate analogous signatures but produce quantitatively different results. Searches for 
cold components in the binary fractions of resonant KBOs are currently practical. The additional 
migration-generated trends described here may be distinguished with objects discovered by LSST. 

Subject headings: celestial mechanics — Kuiper belt: general — planet-disk interactions — planets 
and satellites: formation 



1. INTRODUCTION 



Our Kuiper belt, a ^0.0 1-0. IMm (e.g.JBernstein et al.l 
[200I IChiang et all [20071 : iFuentes fc HolmanI l2nn8n "coh 
lection of planetesimal debris located beyond the orbit 
of Neptune, provides a unique window into the dynam- 
ical processes that shaped the young solar system. The 
orbits of Kuiper belt objects (KBOs) comprise a dynam- 
ical fossil, recording the movements of the giant planets 
during the era of their formation. 

Though almost two decades have passed since the 
discovery of the fir st KBO after Pluto and Charon 
(jJewitt fc Luul Il993f) . and though the orbits of more 
than 500 trans-Ne ptunian objects h ave been well- 
chara cterized fe.g.. lElliot et al.1 l2005t iKavelaars et al.l 
120091 ) ■ a detailed theoretical understanding of the Kuiper 
belt's rich dynamical structure remains elusive. For our 
purposes, we will focus on the following aspects of this 
structure, which remain difficult to reconcile. 

• Approximately 10-20% of all KBOs, including 
Pluto and Charo n, occupy mean mo ti on resonances 
with Neptune (Ij ewitt et al] 119981: iTrujillo et all 
I2OOII : [Kavelaars et al. 2009J), and this fraction may 
be larger if many objects reside in distant reso- 
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nances which are not yet well-characterized. In 
resonance, the orbital period of a KBO forms 
an approximate integer ratio with that of Nep- 
tune, leading to a coherent exchange of energy 
and angular momentum with the planet which is 
periodic over ^10^-10^ year timescales. Though 
reson ance occupat ion enhances long-term stability 
fe.g. jLvkawka fc^ukai 2005i) . the large number of 
resonant KBOs suggests a special dynamical origin 
for this population. 

• Several lines of evidence suggest that classical 
KBOs, those objects that are not in resonance and 
are not curren tly undergoing close e ncounters with 
Neptune (e.g.. lGladman et "aLll2008[ ). consist of two 
populations with different histories and having dif- 
ferent typical inclinations. The inclination distri- 
bution of the belt is bimodal, best fit by a popula- 
tion of bodies with i ^ 5° and a distinct population 
containing equal to greater mass with inclinations 
ranging up to ^35 ° (|Brownll200ll:lElliot et al.ll2005l: 
iGulbis et al.ll201Ci| ). In addition, low and high incli- 
nation classical KBOs have systematically different 
physical properties. The high inclination popula- 
tion is bluer in color (Tcgler fc Rpmanishin 200(1 
iTruiillo fc Brownl l200l iPeixinho et al.ll2008D. has 
a low e r fraction o f wide b inaries (jStephens fc Nolll 
[200l iNoU et al.l l2008bD . and contains all of 
the b rightest (and hence likely largest) KBOs 
([Levison fc Stern ,2001i) . The critical inclination 
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separating those objects with bl uer and redder col- 
ors is ^10° rather than ^^5° (jLvkawka fc Mulffll 
120051: iPeixinho et"aIll2008D . The significance of this 
discrepancy, if any, is not yet known and may be 
complicated by overlap of the low and high inclina- 
tion populations. Throughout this paper, we refer 
to this entire population as the classical belt, the 
low inclination subset as the cold classicals and the 
high inclination subset as the hot classicals. When 
specific numbers are used, we define "low inclina- 
tion" as those objects with i < 10° and "high in- 
clination" as those with i > 10°. We choose 10° as 
our cutoff based on the distributions of colors with 
inclination because we will be investigating a phys- 
ical property of KBOs — their binary fraction. Our 
conclusions will not change substantially if 5° is de- 
termined to be more appropriate for the classical 
beh. 

• Classical KBOs exist with eccentricities ranging up 
to ~0.2. However, cold classical KBOs have sys- 
tematically low eccentricities, likely refiecting a dy- 
namically unexcited population. We note that for 
the purposes of this paper, "cold" and "hot" refer 
only to inclination. 

• The inclinations of KBOs in 3:2 mean motion res- 
onance span the same range as the inclinations 
of the classical belt, but the dynamical data may 
be fit with a sing le high inclination population 
(IGulbis et al.l[20Tnll . 

• At ~48 AU, the densit y of cold classical KBOs 
declines precipit ously (iTruiillo fc BrowiJ 120011 : 
lAUen et al.ll200d) . This transition is located near 
Neptune's 2:1 mean-motion resonance. Though the 
presence of the edge is clear, whether it is an edge in 
semi- major axis corresponding to the 2:1 resonance 
is less certain due to observational selection effects 
at these large distances from the Sun. Further, this 
edge may not be present in the hot classical popu- 
lation (Petit et al., in prep). 

No dynamical model has yet simultaneously reproduced 
the inclination distribution of the belt, its eccentricity 
distribution, and the resonant population of KBOs. 

Two qualitatively different explanations have been pro- 
posed for the abundance of KBOs occupying mean- 
motion resonances with Neptune. The first, which we 
ref er to as "migration-in duced capture," was proposed 
by iMalhotral (|1993l I1995D to explain the orbit of Pluto. 
In this model, KBOs are captured into resonance as Nep- 
tune migrates outward during the late stages of planet 
formation. This migration is powered by planetesimal 
scattering in a d i sk containing ~30Mm in planetesimals 
(jFernandez fc Id 119841: iGomes et ahl I2004D . which also 
acts to damp the eccentricity and inclination of Nep- 
tune. In the simplest version of this scenario, Neptune 
migrates ^10 AU with low eccentricity and inclination, 
though it may well have occupied a more excited orbit at 
earlier times. As this migration progresses, Neptune en- 
trains KBOs in its mean motion resonances as the semi- 
major axes of the reso nances sweep through the primor- 
dial planetesimal disk (IMalhotral 1 19931 11995D . This pri- 
mordial disk must contain some KBOs with pre-excited 



eccentricities (jChiang et al.l 120031 ) to explain the occu- 
pation o f high-order resonance s such as the 5:2. In this 
context, iLevison fc Morbidellil (|200l argue that the en- 
tire Kuiper belt might have been pushed out from closer 
to the Sun by migration, citing the belt's low mass and 
external edge near Neptune's 2:1 resonance. 

As explored bv iHahn fc Malhotral ()2005D . early out- 
ward migration of Neptune reproduces the resonant 
structure of the belt well, but suffers from two diffi- 
culties. First, resonance capture excites the eccentrici- 
ties of captured particles, but does not change their in- 
clinations substantially, so inclinations must be excited 
prior to migration or an alternati ve inclination excita- 
tion mechanism must be found. iGomesI (|2003[ ) argue 
that the large inclination objects in the belt may have 
arisen from scattered objects which interacted with sec- 
ular resonances during migration. This process has a low 
efficiency ^Chiang et al. 20071 . however, making it diffi- 
cult to explain the number of large inclination objects. 
Second, smooth migration does not sufl[iciently excite 
the eccentricities of non-resonant ob jects to explain the 
hot classical population of KBOs. ILevison fc Morbidellil 
(2003) argue that stochastic effects during migration 
could cause KBOs with excited eccentricities to be lost 
from resonance, populating the hot classical belt. How- 
ever, iMurrav-Clav fc Chiang] ()2006l ) demonstrate that 
for a substantial number of KBOs to be lost from res- 
onance during migration, more than a few percent of the 
mass in the planetesimal disk driving migration must be 
in bodies with radii larger than ©(lOOOkm), compara- 
ble to the size of Pluto. This mass fraction is inconsis- 
tent with estimates of the size distributi on in the early 
disk resulting from coagulat ion (e.g., iKenvon fc Luul 
ll999HSchlichting fc Sarill2011l l. Whether this conclusion 
is altered by recent models of planetesimal formation 
tha t may allow rapid formation of 100-1000 km b od- 
ies (jYoudin fc GoodmanI [20051: iJohansen et al.ll2007D re- 
mains to be seen. 

A third difficulty for migration-induced capture arises 
from the different inclination distributions in the classi- 
cal Kuiper belt and the 3:2 resonant population. Res- 
onant capture does not substantially alter inclinations, 
implying that the inclination distribution of a resonant 
population should match that of the population from 
which it captured its members. Though the 3:2 reso- 
nance did not pass through the region of the currently 
observed classical belt, it captured its population from 
an adjacent region of the protoplanetary disk, where one 
would naively expect the inclination distribution to be 
the same. Nevertheless, it is possible that objects caught 
into the 3:2 resonance experienced additional dynamical 
sculpting, rendering this expectation invalid, though no 
model for such sculpting as been produced thus far. In 
particular, we note that if distinct "hot" and "cold" pop- 
ulations exist in the 3:2 resonance, the characteristic in- 
clination dividing these populations may not match that 
in the classical belt. The possibility of additional sculpt- 
ing is less plausible for the 2:1 resonance, which did pass 
through the classical belt. Unfortunately, the inclination 
distribution of the 2:1 is not yet constrained well enough 
to confirm or rule out a cold component. 

A second possibility for the population of mean-motion 
resonances was suggested more recently, inspired by the 
Nice model for the early evolution of the solar system 
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giants ()Tsiganis et al.l [20051 ) . In this model, Jupiter and 
Saturn cross their mutual 2:1 mean motion resonance as 
a result of migration driven by planetesimal scattering. 
This resonance passa ge substantially alter s the orbits of 
Uranus and Neptune. iLevison et alJ (|2008f ) argue that it 
is difficult to place Neptune interior to 21 AU at the end 
of the Nice model, and suggest chaotic population of the 
Kuiper belt as an alternative to migration-induced reso- 
nance capture. They argue that Neptune was scattered, 
with large eccentricity, onto an orbit with a semi-major 
axis close to its current value of 30 AU. Conservation of 
angular momentum during the scattering event leaves the 
planet on a low-inclination orbit. Initially after the scat- 
tering, orbits with semi-major axes ranging from that of 
Neptune out to its 2:1 resonance are chaotic. This chaos 
results from overlapping mean-motion resonances which 
are rendered wide by Neptune's large eccentricity. As 
Neptune's eccentricity is damped by its interactions with 
planetesimals, the resonances decrease in width, ending 
in their current configuration. Objects in mean motion 
resonances are more likely to remain stable during this 
process than their non-resonant counterparts, particu- 
larly at high eccentricities. We refer to this process as 
"chaotic capture." 

We note that this model does have a final phase of 
slow migration of Neptune on a low-eccentricity, low- 
incli nation orbit, howeve r its extent is less than a few 
AU. ILevison et all (pOOl argue that the high inclina- 
tions achieved in their simulations result from this late 
migration phase and the i Gomes (2003) mechanism. The 
inclination distribution of the Kuiper belt is not well- 
matched in this or any model, leaving the source of the 
high inclination population of the Kuiper belt an impor- 
tant outstanding problem. Chaotic capture is appealing 
due to its natural explanation for the edge of the belt 
coincident with Neptune's 2:1 mean-motion resonance. 
It does not, however, reproduce the low-inclination, low- 
eccentricity classical Kuiper belt. Furthermore, an in 
situ belt having the properties of the classical KBOs is 
excited too much by Neptune while the planet's orbit is 
eccentric. 

In this paper, we present a new method for distinguish- 
ing between competing models of Kuiper belt sculpting 
based on observations of the binary fraction of KBOs. A 
significant fraction of large KBOs exist in binaries. More 
than 70 such systems are currently known and this num- 
ber will continue to rise as a result of new KBO searches 
including Pan-STARRS and LSST. INoll et al.l (|2008bf ) 
have shown that the fraction of KBOs that are binaries 
varies with dynamical class in the Ku iper belt. In partic - 
ular, cold classical KBOs (defined in iNoU et al.l ()2008b[ ) 
to have i < 5.5 deg) have a binary fraction of ~ 29%, 
while hot classical KBOs have a fraction of only ~ 10%. 
Furthermore, the physical characteristics of binaries in 
the two pop ulations differ. All of the new binary systems 
reported by iNoll et all ()2008bf ) with low heliocentric in- 
clinations are similar brightness systems and therefore 
presumably consist of roughly equal mass companions. 
On the contrary, among the binary systems in the hot 
population, there were both similar brightness systems 
and systems with lar ge brightness diffe r ences between the 
binary components. iStephens fc Noll! ()2006[ ) found that 
the rate of binaries for the hot classical, resonant, and 
scattered populations combined is about 5.5%. Binary 



statistics are not yet sufficient to establish the binary 
fraction in the resonances alone. 

We hypothesize that this observed variation in binary 
fraction arises ultimately from variations in the location 
in the protosolar nebula at which binary formation took 
place. Broadly speaking, one can identify two classes 
of Kuiper belt binaries. The first class is comprised of 
the large KBOs that are orbited by small satellites. The 
second class consists of similar brightness systems with 
typically wide separations. This first class of systems 
probably originated from a collision and tidal evolution, 
as it has been proposed for t he formation of the Moon 
and the Pluto-Charon system plartmann fc Davi£lll975l : 
iCameron fc War J1976l : lMcKinnonlll989D . Such a forma- 
tion scenario fails however for the second class of Kuiper 
belt binaries, since it cannot account for the large angu- 
lar momenta (i.e. wide separations) of these systems. 

Motivated by this challenge a series of new binary 
formation scenarios, described in fj?! were proposed. 
In many of these scenarios, binaries form by dynam- 
ical capture processes that operated before the Sun's 
planetesimal disk was excited by the giant planets. 
Though dynamical capture and subsequent orbital evo- 
lution can produce binaries with small separations, wide- 
separation, roughly equal mass systems stand out as 
products of dynamical capture rather than collisions, and 
we focus on these systems. In most dynamical capture 
models, the rate at which binaries form varies with he- 
liocentric distance. 

We will demonstrate that migration-induced capture of 
KBOs into mean-motion resonances with Neptune gen- 
crates an observable variation with orbital properties in 
the fraction of KBOs that are wide-separation, roughly 
equal-mass binaries. This variation arises because a cur- 
rent KBO's orbit is correlated with the location in the 
protoplanetary disk at which that KBO formed. In par- 
ticular, we will demonstrate that low-inclination resonant 
Kuiper belt objects could retain a measurable trend in bi- 
nary fraction with heliocentric eccentricity. We calculate 
this trend and argue that if observed, it will provide evi- 
dence both that the cold classical population of the belt 
formed near its current location and that Neptune expe- 
rienced a period of smooth migration. In contrast, if res- 
onant KBOs were emplaced by chaotic capture as Nep- 
tune's eccentricity damped, then resonant KBOs were 
dynamically mixed after excitation of the disk halted bi- 
nary capture, and the binary fraction should not corre- 
late with eccentricity. For narrative clarity, we discuss 
the migration-induced capture mechanism for emplacing 
low inclination resonant KBOs for the majority of this 
paper, and we return to a comparison with chaotic cap- 
ture at the end. This choice should not be interpreted 
as favoring one model above the other. Rather, we hope 
that a future census of Kuiper belt binaries will provide 
a useful discriminant between these classes of models. 

Binary fraction is not the only, or even the first, phys- 
ical property of KBOs that may be searched for dynam- 
ical signatures. In particular, no trend in color with ec- 
centr icity has been ob served for resonant KBOs (e.g., 
[Poressoundiram et al.l (2008). However, for dynamical 
studies, binary fraction has the substantial advantage 
over color that there exists a quantitative theory for how 
it likely varied within the young Sun's planetesimal disk. 
Consequently, we focus on binary fraction throughout 
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the bulk of this paper. We return to a discussion of 
KBO colors and sizes in Section [6] 

In Section [2l we provide a schematic model of 
migration-induced resonance capture. Section [3] de- 
scribes the dynamical history of a KBO captured into 
resonance during a long Neptunian migration. Section |4] 
reviews how binary planetesimals, which will ultimately 
become binary KBOs, form. We calculate formation 
rates as a function of distance from the Sun, and we 
emphasize that this formation must occur before the 
belt is substantially excited or depleted. We combine 
these results to predict the fraction of wide-separation, 
comparable-mass KBO binaries as a function of dynam- 
ical class in Section [5l and we argue that these trends 
should not be present if resonances were populated via 
chaotic capture. Given current projections, these trends 
will be measurable using KBOs discovered by LSST. We 
consider our test in the context of established trends in 
KBO colors and sizes (Section 15]), summarize, and con- 
clude (Section [7]). 

2. FRAMEWORK 

In this section, we assert a schematic of the dynami- 
cal origin of resonant and classical Kuiper belt objects, 
inspired by migration-induced capture. We will use this 
model as a framework for the majority of this paper. 
We emphasize both that this is not a full dynamical 
model and that we are not advocating this scenario as 
the true history of the outer solar system. Rather, we 
will demonstrate that dynamical models which share the 
broad characteristics of our schematic generate a unique 
observable signature in binary KBO dynamics which will 
be absent in other types of models including chaotic cap- 
ture. Thus, it presents a useful framework for interpret- 
ing observations. 

In our scenario, planetesimals interior to ^30 AU from 
the Sun have their eccentricities and inclinations excited 
by the giant planets after the planets themselves undergo 
a dynamical excitation event. The prevale nce of extra- 
solar planets with large eccentricities (e.g.. iButler et afl 
l2006f l suggests that such events are common, and most 
current models of planet formation include dynamical 
upheavals, whether t hey result from reso nance crossings 
as in the Nice model ()Tsiganis et al.|[2005l ) or from reduc- 
tions in the efficiency of damping by small bodies and gas 
(jGoldreich et al.l 12004a'). Beyond --30 AU, a planetesi- 
mal disk formed and survived with low eccentricity and 
inclination. After this period of dynamical excitation, 
the eccentricity and inclination of Neptune were damped 
by the residual planetesimal disk at a semi-major axis of 
~20-25 AU. It then migrated outward to 30 AU, main- 
taining low eccentricity and inclination, as it scattered 
planetesimal debris (jFernandez fc IplllQSl ). 

Whether the above scenario can reproduce, in detail, 
the eccentricity and inclination properties of the belt re- 
mains to be determined and is the subject of ongoing 
research. For our purposes, the relevant characteristics 
of a model such as this are as follows. The cold classical 
belt, with low inclinations, formed in situ. High inclina- 
tion KBOs formed closer to the Sun than low inclination 
KBOs. The high inclination and low inclination pop- 
ulations of the Kuiper belt were generated by distinct 
dynamical mechanisms, have distinct histories, and were 
formed at different locations in the disk. Furthermore, 



high inclination KBOs experienced more dynamical mix- 
ing than low inclination KBOs. As Neptune migrated out 
through the planetesimal disk, it captured both high and 
low inclination objects into resonance. Low inclination 
objects were likely to be captured near their formation 
locations, while high inclination objects may have been 
moved substantially before capture. 

Even in smooth migration scenarios, the hot popula- 
tion of classical KBOs is likely to be dynamically mixed 
from its formation location and may be subject to en- 
hanced binary disruption due to scattering by Neptune. 
Given these considerations, we suggest that the searches 
for coherent signatures of migration in the physical prop- 
erties of KBOs should be restricted to the low inclination 
population. 

3. ORBITAL PROPERTIES OF RESONANT KBOS 
CAPTURED VIA MIGRATION 

The final eccentricity Cfin of an object caught into p:q 
resonance by a slowly migrating Neptune can be related 
to the semi- major axis at which it was captured, Ccap, 
using an adiabatic invariant of the system. After review- 
ing relevant properties of resonances, we describe this 
invariant, originally derived for a system containing a 
single planet on a circular orbit, interacting with a test 
particle occupying the same plane (Section 13. demon- 
strate that it is useful under a range of likely histories 
for the outer solar system (Section 13. 2p , and comment 
on the relationship between eccentricity and inclination 
among observed resonant KBOs (Section 13. 3|) . 

3.1. Orbital Preliminaries and Brouwer's Constant 

First, we briefly review pertinent facts about mean mo- 
tion resonances. Further discussion may be found in, 
e.g., the textbook Murray & Dcrmott (2000). Consider 
a system consisting only of a central star with mass M, , 
a planet with mass Mp on a circular orbit, and a col- 
lection of test particles. A particle in external p:q mean 
motion resonance with the planet experiences libration 
of the resonance angle 

(j> = pX-qXp- {p- q)zD , (1) 

where A and Ap are the mean longitudes of the test par- 
ticle and the planet, respectively, and w is the longitude 
of pericenter of the particle. This expression is appro- 
priate for resonances dominated by the eccentricity of 
the particle, as is the case for all known resonant KBOs. 
For a first-order resonance {p — q = 1), the angle (f) is ap- 
proximately the angular position of the test particle with 
respect to pericenter at the time of conjunction, when the 
planet lies on the line connecting the particle to the star. 
Because the gravity of the planet generates the largest 
perturbation to the particle's orbit at conjunction, if (f) 
remains constant, then the planetary perturbations co- 
herently add. This resonant effect can generate large 
changes to the particle's orbit. 

When resonance is not exact, (f) is not a con- 
stant, but rather librates through a limited range of 
values. The tiniescale for this libration is inb ~ 
Ciiorb[(M*/Mp)e-|P-«l]i/2^ where the particle has or- 
bital period torb = 27r/f2, orbital angular frequency fl, 
and eccentricity e. Libration within the resonance pro- 
duces variations in the particle's semi-major axis, a, 
with a maximum magnitude (full width) of max(i5aiib) ~ 
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2C2a[(Mp/M,)elP-'?l]i/2. The coefficients Ci and C2 are 
order unity constants. For the 3:2 resonance, Ci ~ 0.3 
and C2 ~ 3.6. For the 2:1, Ci « 0.9 and C2 ~ 3.0. 

Now imagine that the planet's orbit expands, while 
remaining circular. As long as this outward migration is 
sufficiently slow, the orbit of a particle in external mean 
motion resonance also expands, maintaining the resonant 
lock. As it does so, its eccentricity evolves such that 



Cb - a{VT 



(2) 



where Cb, known as Bro uwer's c onstant, is an adiabatic 
invar i ant of the system (IBrouwe r 1963; Yu & Trcmainc 
119991: IHahn fc Malhotrall2005i r^his relationship applies 
as long as the particle remains in resonance and the 
timescale of migration tniig = flp/op ^ tsyn, the syn- 
odic period (time between conjunctions). In general, the 
former is a stricter constraint. To capture and maintain 
particles in resonance, the planet must migrate slowly 
enough that [max((5aiib)/a]tinig ^ f^b. 

Fundamentally, Brouwer's constant encodes the ratio 
between angular momentum and energy transfer during 
transport in resonance. For a resonant particle, the spe- 
cific energy E = — GM,/(2a) and angular momentum, 
L ~ (GM,a(l — e^))^/^, measured with respect to the 
host star, increase at rates E and L, respectively, such 
that 

E w Lrip (3) 



when averaged over a synodic period fe.g.. iChiang et all 
I2007I) . In contrast, maintenance of a circular orbit would 
require E = Lfl. Equation ([3]) may be intuitively un- 
derstood as the result of forcing at the planet's orbital 
frequency flp. The radius f and the velocity w of a Ke- 
plerian orbit are roughly perpendicular to one another, 
and their magnitudes differ by approximately fl, so one 
might expect a perturbing force F to generate a torque 
•rxF and rate of work F ■ v that differ by approximately 
This intuition fails because the majority of the energy 
and angular momentum change averages to zero over a 
synodic period. Only second-order effects remain, which 
instead yield a term proportional to the forcing frequency 

In this paper, we appeal to Brouwer's constant to relate 
the semi-major axis at which a particle was caught into 
resonance during Neptune's migration to its current ec- 
centricity in resonance. For example, a straightforward 
application of Equation ([2]) implies that under the hy- 
pothesis of smooth migration, a KBO currently having 
Cfin = 0.2 and semi-major axis = 39.5 AU in Nep- 



tune's 3:2 resonance was captured at Oc 



34.9 AU if 



it initially had zero eccentricity. This consideration is 
the basis for estimates of the extent of Neptune's pri- 
mordial migration — if KBOs were all caught with low 
eccentricities, the current eccentricities of resonant ob- 
jects imply that Neptune began its migration at least 10 
AU closer to the Sun than it s present orbit (Malhotra 
[1995 IHahn fc Malhotol l2005l) . In Section E^l we dis- 
cuss the accuracy of this mapping in the context of the 
solar system. 

3.2. Conservation of Brouwer's Constant Under 
Realistic Migration Conditions 



Given the assumption that smooth migration of a low- 
eccentricity Neptune occurred, one may still worry that 
difficulties will arise in applying Brouwer's constant, de- 
veloped for the circular restricted three-body problem, to 
the true solar system. The signal might be washed out 
due to the presence of four giant planets, all with small 
but notable eccentricities and inclinations, by the fact 
that we do not, a priori, know the initial eccentricities 
of KBOs, or by chaotic evolution over the full age of the 
solar system. Fortunately, this is not the case. 

To evaluate our ability to use Brouwer's constant to 
identify the semi-major axis at which a KBO was cap- 
tured into resonance, we performed a numerical integra- 
tion of KBOs evolving under the influence of migrating 
planets. We performed this integration in two steps, both 
using the hybrid pa ckage of the N-b ody integrator Mer- 
cury, version 6.2 ( Cha mbers! Il999f l. In the first step, 
we begin with 10,000 test particles with initial eccen- 
tricities ranging uniformly between and 0.2 and initial 
semi-major axes ranging uniformly between 25 and 50 
AU. We include all four giant planets with their current 
eccentricities and inclinations. The initial semi-major 
axes of the planets differ from their current values by 
Aa = 7, 3, 0.8, and -0.2 AU for Neptune, Uranus, Sat- 
urn, and Jupiter, respectively, and they evolve to their 
current semi- major axes over an exponential timescale 
r = 10^ as in iMalhotral ()1995[ ). This migration is imple- 
mented by applying accelerations to the giant planets of 
the form 0.5(d/a)v in the Mercury6 routine mfo_user, 
where d/a = (Aa/a)T~^ exp(— </t) at time t and v is 
the vector velocity of the particle. This integration lasts 
for 3 X 10^ years with a time step of 8 days. At the 
end of 3 X 10^ years, we continue the integration with 
no migration force and with a time step of 150 days for 
a further 10^ years only for the 1350 objects that were 
captured into 3:2 (647) or 2:1 (703) resonance with Nep- 
tune. Because we do not employ the migration force over 
the entire integration, the planets end with semi-major 
axes that differ from their true values by approximately 
exp(— 3)Aa w 0.05Aa. After 10^ years, 65 objects re- 
main in the 3:2 resonance and 47 in the 2:1. We note that 
this corresponds to a retention efficiency smaller than the 
39 % and 24% retention efficienci es calculated over 1 Gyr 
bv lTiscareno fc Malhotral (j2009( ) for the 3:2 and 2:1 res- 
onances, respectively, due to chaotic diffusion under the 
infiuence of the four giant planets. This difference likely 
results either from differences in the population of reso- 
nant objects at the beginning of our long term integra- 
tions or from the slightly different configuration of the 
giant planets in our long-term integration. Our integra- 
tions are intended to verify our ability to retrieve initial 
semi-major axes under a wide range of initial conditions 
rather than to faithfully reproduce the initial conditions 
in the disk, so we do not consider our substantially lower 
retention fraction significant. 

Figure [1] compares the true initial semi-major axis with 
the initial semi-major axis Ocap calculated using Equation 
^ for the KBOs which were captured and maintained in 
3:2 or 2:1 resonance for the full 10^ years. In order to ap- 
ply Equation ([2]) without a priori knowledge of the initial 
orbits of the particles, we assume zero initial eccentricity 
so that 

afin rr^ ■ W 
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Figure [T] exhibits a clear trend between the true initial 
a and the calculated value. The scatter in this relation 
corresponds to a cut at e = 0.2 misclassifying the ini- 
tial semi-major axis of approximately 20% of observed 
objects. This scatter can be substantially reduced by 
evaluation of a second adiabatic invariant of the system 
which permits us to estimate the initial eccentricity of 
a captured KBO, allowing for a more accurate calcula- 
tion of the initial semi-major axis of a captured particle 
(Murray-Clay, Lithwick, & Chiang, in prep.). However, 
this improved calculation is substantially more compli- 
cated than application of Equation ([2]) and is beyond 
the scope of this paper. We have verified that our abil- 
ity to recover the initial semi-major axis is not lost after 
evolution over the age of the solar system. 
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Fig. 1. — True initial semi-major axis of 3:2 (plusses) and 2:1 
(circles) KBOs captured into resonance via migration in a numer- 
ical simulation compared against the semi-major axis at the time 
of capture inferred from Equation l(4]l. The integration includes 
all four giant planets with their current eccentricities and inclina- 
tions. Test particles have initial eccentricities ranging from to 
0.2. This trend may be substantially tightened by consideration of 
a second adiabatic invariant of the system (Murray-Clay, Lithwick, 
& Chiang, in prep). 

Since all particles in a given resonance have compara- 
ble semi-major axes. Equation (j?]) implies that particles 
caught into resonance at smaller semi-major axes have 
larger final eccentricities and that resonant evolution can 
lead to substantial eccentricity excitation. This excita- 
tion was one of the o riginal motivati ons for the idea of 
Neptune's migration (IMalhotralll993h . 

3.3. Eccentricity vs. Inclination 

We note that differences in the relationship between 
eccentricity and inclination, i, for non- resonant and res- 
onant objects are broadly consistent with migration- 
induced capture. Figure [2] displays the relationship be- 
tween e and i for known classical and 3:2 resonant KBOs 
in the Minor Planet Center (MFC) catalog with at least 



3 observed oppositions. Inclinations are calculated with 
respect to the invariable plane of the solar system. The 
invariable plane differs from the true dynamical plane 
of the Kuiper belt — known as th e forced or Laplace 
plane — by at most a few degrees (jBrown fc Panl[200l 
.Chiang fc Choill2008f). We use the dynamical classifica- 
tions of iGladman et al.l ()2008[ ). which includes objects 
with 3 oppositions as of May 2006. For the approxi- 
mately 20% of each population of objects not contained 
in that work, we assign a tentative dynamical classifica- 
tion based on a 10^ year integration of the nominal orbit 
in the presence of the four giant planets0 
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Fig. 2. — Eccentricity vs. inclination with respect to the invari- 
able plane for observed classical KBOs (black crosses) and 3:2 reso- 
nant KBOs (red circles) from the MFC catalog. The eccentricities 
of 3:2 resonant KBOs are systematically higher than the eccen- 
tricities of classical KBOs and do not correlate with inclination. 
This pattern is consistent with migration-induced eccentricity ex- 
citation, which does not affect inclination. Because 3:2 objects are 
more stable than classical bodies at high eccentricity, this pattern 
is also consistent with removal of high e classical objects due to 
planet scattering, as occurs, for example, in chaotic capture. We 
include for reference the line e = \/2sin(i), representing equiparti- 
tion excitation of eccentricity and inclination. 

As mentioned in Section [TJ 3:2 objects do not exhibit 
the same excess of low inclination bodies seen in the clas- 
sical population. This lack of a large cold population 
is potentially problematic in the context of migration- 
induced capture, but is not conclusive because the 3:2 
resonance did not sweep through the region of the clas- 
sical belt. However, among hot objects, the distribution 
of eccentricity and inclination mimics that expected for 
migration-induced capture. The inclinations of KBOs in 
the 3:2 resonance are not correlated with eccentricity and 
span the same range as the inclinations of hot classical 
objects. At the same time, the eccentricities of the res- 
onant KBOs are systematically higher than those of the 
classical objects. This pattern of similar inclinations but 

^ Objects for which the 3:2 resonance angle <f> librates are as- 
signed to the 3:2 resonance. An object is designated classical if 
its semi-major axis a > 30 AU, its eccentricity e < 0.24, it sur- 
vives the 10^ year integration, <j> does not librate for the 5:4, 4:3, 
3:2, 5:3, 7:4, 9:5, 2:1, 7:3, or 5:2 resonances, and finally either 
its pericenter q = a{l — e) > 35 AU or the Tisserand parameter 

T = a iv / d -|- 2 \/ (a/ajv)(l — e^) cos(Ai) > 3, where ajv is the semi- 
major axis of Neptune and Ai is the mutual inclination between 
the planes of Neptune and the KBO. 
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offset eccentricities is consistent witti migration-induced 
capture, wliicli excites eccentricities wfiile leaving incli- 
nations largely unchanged. This pattern is suggestive 
but not conclusive evidence for migration-induced cap- 
ture because high eccentricity KBOs are more stable in 
the 3:2 resonance than in the classical belt, and the ob- 
served eccentricity offset could alternatively result from 
the removal of high e classical bodies. 

4. BINARY FORMATION AS A FUNCTION OF 
SEMI-MAJOR AXIS 

Having established that under the migration-induced 
capture scenario, we can statistically recover the the ini- 
tial semi-major axes of KBOs in mean motion resonance 
with Neptune, we now show that the binary fractions 
of different populations of KBOs are expected to be a 
function of the locations in the protoplanetary disk at 
which those populations formed. In Section 14.11 we re- 
view proposed binary formation scenarios. We argue that 
binary formation mediated by dynamical fricti o n wit h 
small bodies, as proposed by iGoldreich et al.l ()2002l), 
likely formed the Kuiper belt binari es. In Section 14.21 
we review the IGoldreich et al.l (j2002| ) formation mecha- 
nism, and in Section [4.3| we calculate the binary fraction 
as a function of formation distance. We note that other 
binary formation mechanisms which vary with distance 
could also generate an observable, though different, sig- 
nature in the resonant population. 

4.1. Binary Formation Scenarios 

High mass-ratio binaries in the Kuiper b elt, includ- 
ing Pluto/Charon, likely formed in collisions (iMcKinnonI 
[198 9). However, the majority of observed binary 
KBOs are roug hly equal-mass, wide-separation binaries 
(|Noll et al.ll2008M that have too much angular momen- 
tum to be formed by the same mechanism. Motivated by 
this challenge, a series of new binary formation scenarios 
was proposed (e.g. fWcidenschilling 2002; Goldrcic h et al. 
2002; Funato et al. 2004; Astakh ov et al.ll2005l:lLee et al. 
2007; Nesvorny ct al. 2010). Most of these scenarios ap- 
peal to interactions within the Hill sphere, the region 
interior to the Hill radius of a KBO. The Hill radius de- 
notes the distance from a KBO at which the tidal forces 
from the Sun and the gravitational force from the KBO, 
both acting on a test particle, are in equilibrium. It is 
given by 

where a is the semi-major ax is of the KBO and m its 
mass. IWeidenschillind ()2002D proposed a colhsion be- 
tween two KBOs inside the Hill sphere of a third. How- 
ever, in the Kuiper belt, gravitational scattering be- 
tween the two intruders is about 100 time^ more com- 
mon than a collision. Binary formation by three body 
gravitational de fl ection (L^ mechanism), as proposed by 
IGoldreich et al.l (|2002l ). should therefore dominate over 
such a collisional formation sce nario. A second binary 
formation scenario proposed by IGoldreich et al.l ()2002l ) 
consists of the formation of a transient binary, which 

® For this estimate we used that a ~ R/Rjj ~ 10^'' and assumed 
that the velocity dispersion of the KBOs at the time of binary 
formation is less than their Hill velocity (see il4.2|l . 



becomes bound with the aid of dynamical friction from 
the s urrounding sea of s mall b odies (L^ mechanism, 
S3)- iSchhchting fc Saril ()2008al ) demonstrated that at 
the distance of the K uiper belt, domiii ates over L^. 
'Asta khov et al.l (|2005H and lLee et alj (|2007[ ) suggest that 
transient binaries that spend a long time in their mutual 
Hill sphere, near a periodic orbit, form the binaries in 
the and mechanisms. iSchlichting fc Sarll ()2008al ) 
investigated the relative importance of these long-lived 
transient binaries. They found that such transient bi- 
naries are not important for binary formation via the 
mechanism and that they become important only 
for very weak dyna mical friction in the L^s mechanism. 
iFunato et al.l ()2004D proposed a binary formation mech- 
anism which involves a collision between two large KBOs 
that creates a small moon. An exchange reaction replaces 
the moon with a massive body with high eccentricity and 
large semi - major axis. 

Finally, INesvornv et al.l (j2010t) suggested that Kuiper 
belt binaries could form directly during a gravitational 
collapse that leads to the formation of 100-km-sized 
KBOs. This binary formation scenario is fundamentally 
different in the sense that it assu mes, based on recent 
work on the streamin g inst ability (|Youdin fc GoodmanI 
I2005t IJohansen et"an i2007D . that 100-km sized KBOs 
formed by direct gravitational collapse rather than coag- 
ulation. Whether such a gravitational collapse KBO for- 
mation scenario is feasible and whether it can explain the 
observed KBO size distrib ution, which is w^ell matched by 
coagulation simulations (iKenv on fc Luul Il999t iKenvonI 
120021 : IKenvon fc Bromle v 2004), remains to be deter- 
mined. Furthermore, a gravitational collapse binary for- 
mation scenario, would have to provide a convincing ex- 
planation for the large binary fraction of similar bright- 
ness systems in the Kuiper b elt and the absence of such 
binaries in the Asteroid belt. INesvornv et al.l (|2010[ ) ap- 
peal to enhanced binary disruption in the Asteroid belt 
due to collisions and scattering events to explain this dif- 
ference. A dynamical origin, which is common in all the 
other binary formation scenarios discussed above, takes 
advantage of the fact that the Hill radius is more than an 
order of magnitude larger for KBOs compared to similar 
sized objects in the Asteroid belt. Gravitational capture 
binary formation scenarios therfore support the obser- 
vations that the Kuiper belt is the ideal place for wide, 
similar brightness binaries to form. 

We decided to focus on the IGoldreich et al.l ()2002D bi- 
nary formation mechanism for the calculations in this 
paper. Our choice is motivated by the success this bi- 
nary formation scenario has had in explaining the bi- 
nary abundance and the observed binary separation dis- 
tribution. It also predicts the existence of triplet sys- 
te ms. The first t riplet system was reported recently 
bv iBenecchi et ahl (|2010[ ). One set of o bservations that 
are n ot successfully explained by the IGoldreich et al.l 
(|2002| ) formation scena rio is the mu tual binary inclina- 
tion distribution (Noll et al.l l2008al ) . Binary formation 
rates are extremely sensitive to the velocity dispersion 
V ~ (e) ria of the population of large planetesimals from 
which they form, where (e) is the typical eccentricity 
of a planetesimal. Efficient binary formation requires 
that the large KBOs have a velocity dispersion that is 
less than their Hill velocity, vh = flRn- Binary for- 
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mation rates quickly exceed the age of the solar system, 
once the velocity dispersion grows significan tly above vh 
(|Noll et al.ll2008al:ISchlichting fc Sarili2008al) . Binary for- 
mation while sub-Hill velocities prev ailed, however, pre- 
dicts low mutual binary inclinations fCh iang et al.|[2007l : 
[Schlichti ng fc Sari 2008b, ) and is inconsistent with the 
rou ghly uniforrn inclin ation distribution that is observed 
(see lNaoz et alj ()2010[ ) and references theirin) . Evolution 
of the mutual binary inclination after formation provides 
a possible solution to this problem. Although it has, 
for example, been suggested that the Kozai mechanism 
coul d affect the orbital e volution of Kuiper belt bina- 
ries (jPerets fc Naod [20091 ) . it remains a topic of ongoing 
research if mutual binary inclinations can be excited to 
sufficiently large values that would allow such a mecha- 
nism to operate. 

4.2. Binary Capture Mediated by Dynamical Friction 

After specifying our assumptions, we now review the 
binary formation rate for the L^s mechanism, which 
will be useful in calculating the exp ected binary frac- 
tion. Following iGoldreich et all (12002) . we use the "two- 
group approximation" , which consists of the identifica- 
tion of two groups of objects, small ones, that contain 
most of the total mass with surface mass density a, and 
large ones, that contain only a small fraction of the to- 
tal mass with surface mass density Y, <^ a. We as- 
sume that a ~ (740 (a/40 AU)^ which is a generaliza- 
tion of the minimum-mass solar nebula (jWeidenschillinel 
[19771 IHavashi|[T98l . where (740 ~ 0.3 g cm ^ is the sur- 
face density at a heliocentric distance of 40AU. The 
power-law index 7 is typically assumed to have values 
between -1 and -1.5, with submillimeter observations of 
the outer regions of pr otoplanetary disks favoring 7 « — 1 
([Andrews et all 120091 ). 

We likewise parameterize the mass surface density 
of large bodies with sizes of i? ~ 100 km and larger 
as E ~ I]4o(a/40 AU)^ , where we treat /3 as a free 
parameter. Estimates from Kuiper belt surveys (e.g., 
Truiillo et al.|[200l iTruiillo fc Brow^[2003[: [Petit et alj 



2008bl: iFraser fc KavelaarsI I2009t iFuentes et all I2009I) 
y ield E4n ^ 3 x 10~^ g c m~^ in the current belt. 

[Goldreich et al.l ([2002D demonstrate that, under the as- 
sumption that E was the same as it is now during the 
formation of Kuiper belt binaries, the velocity dispersion 
V of ~100km KBOs is damped by dynamical friction from 
the sea of small bodies such that v < vh- This is referred 
to as the "shear-dominated velocity regime" because un- 
der these conditions, the relative velocity of two large 
bodies that encounter one another is dominated by their 
Kepleria n shear. Binary formation is inefficient wh en 
v>VH ([Chiang et al.l[2007[: [Schhchting fc Sarill2008bl) . 

In particular, in this scenario, large bodies grow by 
the accretion of small bodies. Large KBOs viscously stir 
the small bodies. This stirring increases the small bodies' 
velocity dispersion, u, which grows on the same timescale 
as R and is given by 



(6) 



^ ^ ^ Uo Au) 

where a = R/Rh (IGoldreich et al.[ [2001 11 Meanwhile, 



the velocity v of large KBOs increases due to mutual vis- 
cous stirring, but is damped by dynamical friction from 
the sea of small bodies such that v < u. Balancing the 
stirring and damping rates of v and substituting for u 
from Equation [51 we find 



3(/3-7)-l-2 
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A transient binary forms when two large KBOs pene- 
trate each other's Hill sphere. This transient binary must 
lose energy in order to become gravitationally bound. In 
the mechanism, the excess energy is carried away by 
dynamical friction with the sea of small bodies. Since 
it is not feasible to examine the interactions with each 
small body individually, their net effect is modeled by 
an averaged force which acts to damp the large KBOs' 
non-circular velocity. We parameterize the strength of 
this damping by a dimensionless quantity D defined as 
the fractional decrease in non-circular velocity due to dy- 
namical friction over time rt~^: 



pR \VH 
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where p 1 g/cm^ is the internal density of a KBO. 
The first expression is an estimate of dynamical friction 
by the sea of small bodies for u > vh- The second ex- 
pression describes the mutual excitation among the large 
KBOs for V < Vh- These two expressions can be equated 
since the stirring among the large KBOs is balanced by 
the damping due to dynamical friction. Using this pa- 
rameterization for the dynamical friction, the binary for- 
mation rate for equal mass bodies via the L^s mechanism 
in the shear-dominated velocity regime is given by 



FR 



L2, 



= At,2,D 



(9) 



([Goldreich et al.ll2002fl where At2, is a cons tant with a 
value of about 1.4 ([Schlichting fc Sariir2008aD . FRl2, ~ 



1.3 X 10~^yr~^ when evaluated at 40 AU. 

4.3. Binary abundances 

We now calculate the expected variation in the fraction 
of binaries formed as a function of semi-major axis. The 
binary separation, s, shrinks due to dynamical friction 
at a rate s~^\ds/dt\. The probability of finding a large 
KBO in a binary per logarithmic band in s is given by 



s p{s) 



FR 



s~'^\ds/dt\ 



(10) 



Collisions among the small bodies and the associated coUisional 



where p{s) is the probability density of finding a large 
KBO in a binary with binary semi-major axis s. There 
are two regimes in binary separation that need to be con- 
sidered separately. In the first regime the binary semi- 
major axis is sufficiently large such that the velocity of 
the binary components around their common center of 
mass, vb, is small compared to the velocity dispersion 
of the small bodies that provide the dynamical friction, 
i.e. vb < u. In this case s~^\ds/dt\ = D, given by 

damping of their velocity dispersion is unimportant on the viscous 
stirring timescale and binary formation timescale, provided that 
the small bodies have radii of order 100m or larger. 



Equation (|S]). However this regime only applies to bi- 
nary semi-major axes that satisly GM/v? < s < Rh 
which translates into binary angular separations of 3" 
and greater at heliocentric distances of about 40 AU. To 
our knowledge, only one (2001 QW322. |Petit et al.|[2009l) 
out of about 70 currently known KBO binaries satisfies 
this criterion. We therefore focus from here onwards on 
the second regime for which vb > u, which applies to 
binaries with angular separation of less than about 3" at 
a distances of about 40 AU. In this regime, the binary 
separation shrinks due to dynamical friction at a rate 



a I u 

pR \VH 
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Substituting the above expression and the binary forma- 
tion rate from Equation ([9|) into Equation (fTO| we find 
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where in the last step, we used the expr ession for u given 
in Equation ^ (jGoldreich et all 120021 ). The only semi- 
major axis dependence that enters in Equation re- 
sults from <T and a. For cr cx the binary fraction as a 
function of a is given by 



s p{s) 



pR 
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at the time of formation. 

Equations ([T2l) and ([T3| show that the binary fraction 
does not explicitely depend on E. We also note that, 
because we are interested in the scaling of s p{s) with 
a, the values of 1740 and S40 do not matter as long as 
their ratio is such that binary formation can proceed in 
the shear dominated velocity regime. Given everything 
else equal, we find that the binary fraction is between 1.4 
times (7 — —3/2) to twice (7 — —1) as high at 50 AU 
compared to 25 AU during the epoch of binary forma- 
tion. This trend in the binary fraction as a function of 
heliocentric distance is in agreement with observations 
that find a binary fraction that is more than a factor of 2 
lower in the excited and hot classical populations, which 
most likely formed closer to the Sun and were scattered 
to their current locations, than in the cold classical pop- 
ulation, which most likely formed in-situ at about 40 AU. 
An analogous calculation for the L'^ mechanism of binary 
formation yields the scaling relation s p{s) oc a'^~^^^, 
where we used E/cr - a3/4 (ISchlichting fc Saril[20lTI) . 

In Section[51 we explore the expected binary fraction as 
a function of current dynamical class in more detail. We 
recall here that the power-law index of the small bodies' 
mass surface density, 7, is uncertain. Often 7 = — 1.5 is 
used based on the minimum mass solar nebula (jHavashil 
|1981|) , but protoplanetary disk observations at dist ances 
comparable to the Kuiper belt ([Andrews et al.|[2009l ) sug- 
gest 7 = — 1. In the calculations that follow, we assume 
7=— 1. If 7 = —1.5 then the trends in the binary abun- 
dance will still be present but sm aller in magnitude. Par- 
ticle pile-ups fe.g.- lYoudin fc Chia ng 2004), if they occur, 
result in a shallower mass surface density profile, making 
trends in binary abundance more pronounced. 

Finally, as mentioned above, efficient binary formation 
(as proposed bv .Weidenschilling, (,2002. ): .Goldrcich et al.l 



(|200^: iFunato et al.l (|200l : lAstakhov et"all (|2005l) : 
ILee et al.l (120071) 1 requires v < vh, which implies that 
the velocity dispersion of the large bodies that form 
binaries has to be damped to sub Hill velocities. If 
this damping is provided by dynamical friction gener- 
ated by small bodies, as we have assumed in deriving 
Equation ((T]), than this implies that S ^ cr, because 
V > Vh otherwise. These binary formation scenarios 
are therefore inconsistent with reduction of the mass in 
the Kuiper-belt by two orders of magnitude or more in 
an entirely size-independent fashion, for example in one 
large scattering event. Analytic work and numerical co- 
agulation simulations that model the growth of KBOs 
indeed find that only about a'^/^ ^ 10^'^ of the initial 
disk mass is converted into large KBOs suggesting that 
E/g y 10^ dKeny on & Luu 1999; Kenvon & Bromle^ 
120041: ISchlichting fc Saril 1201 ID . which is consistent with 
the surface density values used in this section and the 
conclusion that v < vr- 

5. KBO BINARY FRACTIONS AS A FUNCTION OF 
DYNAMICAL CLASS 

In order to test the nature of Neptune's migration 
with the abundance of comparable-mass, wide-separation 
KBO binaries, we suggest the following procedure. We 
recall that we are searching for a signature of an outer 
solar system history in which the cold (low inclination) 
classical population formed in situ, the hot (high incli- 
nation population) formed closer to the Sun and was dy- 
namically mixed during transport to its current location, 
and resonant KBOs are a superposition of these two pop- 
ulations, caught into resonance during migration of Nep- 
tune spanning between several and ~10 AU {^^. Our 
proposed signature of this class of migration histories is 
summarized in Figure [3] 

5.1. Looking for a Cold Component 

First one needs to establish whether a low inclina- 
tion component, analogous to the cold classical pop- 
ulation, exists in a given resonant population. Re- 
call that because any hot population will contain low- 
inclination bodies, a cold component refers not merely 
to objects with low-inclinations but rather to a separate 
low- inclination population. As yet, for the 3:2 there is 
no dynamical evidence for two inclination components, 
implying that if one exists it is less pronounced than in 
the classical belt, and insufficient data is av ailable for the 
2:1 r esonance to make this determination ()Gulbis et al.l 
I2OIOD . We note that because 3:2 resonant KBOs orig- 
inated in a different location in the disk than classical 
KBOs if migration-induced capture occurred, and be- 
cause the inclinations of resonant bodies evolve more over 
long-term integrations than their classical counterparts 
(Volk & Malhotra, in prep.), the characteristic inclina- 
tion that divides the low and high inclination popula- 
tions in the 3:2 resonance could be different from that 
in the classical belt. Furthermore, this characteristic in- 
clination could vary between resonances. If a resonance 
does contain a cold component with a broader width or 
a smaller relative contribution than seen in the classical 
belt, it should be observable in binaries, just as the cold 
component can be picke d out just from t he binary frac- 
tion in the classical belt ()Noll et al.ll2008bD . We therefore 
suggest searching for a cold component by comparing the 
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binary fractions of low and high inchnation objects in a 
given resonance. We discuss this distinction in the con- 
text of KBO colors and sizes in Section [6l 

TheoreticaUy, if cold and hot components are present 
in a resonance, our schematic migration scenario predicts 
fewer binaries in the high inclination population. This is 
both because they originated closer to the Sun and be- 
cause the excitation process may have broken b i naries 
that were originally there. iParker fc Kavelaari ()2010f ) 
have shown numerically that binary disruption is com- 
mon for objects transported across the semi-major axis 
of Neptune. They consi der the particu l ar ch aotic cap- 
ture scenario modeled in iLevison et al.l (|2008[ ) in which 
all KBOs, including those ultimately on low eccentric- 
ity and low inclination orbits, begin interior to the cur- 
rent cold classical belt, and they find that wide binaries 
are efhciently destroyed. This finding is in confiict with 
the o bserved large bin ar y fraction in the cold clas sical 
belt (jNoll et al.ll2008al ). IParker fc KavelaarsI (|20Tol) ar- 
gue that the difference in binary fraction between the 
cold and hot classicals may reflect scattering of the hot 
classicals off of Neptune during their history, contrasted 
with a gentler dynamical history for low inclination ob- 
jects, which formed exterior to Neptune and never en- 
countered the planet. Depending on the details of the as 
yet unknown process exciting inclinations in the Kuiper 
belt, there may be trends with inclination in the binary 
fraction of the high inclination population. Nevertheless, 
we expect the overall fraction to be smaller than that in 
the low inclination population. 

A lack of a cold component in the 3:2 resonance might 
be explained in the context of migration-induced cap- 
ture if the original disk was substantially excited interior 
to the current location of the 3:2, though no proposed 
model currently produces such excitation. In contrast, 
a lack of a cold component in the 2:1 resonance would 
seriously challenge a scenario invoking smooth migration 
through a hot plus cold planetesimal disk because the 2:1 
resonance should have swept through the location of the 
current cold classical belt. 

We emphasize that the 3:2 and 2:1 resonances are the 
best locations to test for the presence of a cold com- 
ponent. Higher order resonances preferentially capture 
objects with large eccentricities. Therefore, if eccentric- 
ities and inclinations are correlated, higher order reso- 
nances may only contain hot components even if they 
swept through a disk that had a cold component. How- 
ever, if a higher order resonance does have a population 
of low inclination objects, it should also have a higher 
binary fraction in the cold component compared to the 
hot component. 

5.2. Looking for Trends in the Low Inclination 
Resonant Population 

Provided that a cold component can be established in 
a Kuiper belt resonance fSection [5?T|) . we suggest looking 
for the following trends in binary fraction within the low 
inclination resonant population. First, we propose look- 
ing for differences across the binary fractions of the cold 
populations of different resonances and of the cold clas- 
sical belt, and second we propose looking for trends in 
binary fraction with eccentricity within the low inclina- 
tion population of a single resonance. Because it requires 
observations for fewer objects per resonance, the former 
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Fig. 3. — Fraction of low-inclination, resonant Kuiper belt objects 
predicted to be wide-separation, roughly equal-mass binaries, un- 
der the hypothesis that Neptune captured objects into resonance 
as it migrated ~10 AU. We normalize this expected fraction to 
the binary fraction among low-inclination (cold) classical objects 
(CKBOs) at 45 AU. The vertical lines show both the maximum 
eccentricity excitation from migration by 9.4 AU (calculated using 
Equation and the maximum observed eccentricities in 3:2 and 
2:1 resonances. The agreement of these values is promising for the 
migration scenario, but chaos at large eccentricities might also ex- 
plain the different maximum eccentricities in different resonances. 
The observed binary fraction among highly inclined (hot) CKBOs 
(shaded region) is consistent with formation <15 AU from the Sun, 
though scattering events could disrupt binaries, complicating this 
interpretation. 



effect may be easier to measure, and we describe it first. 

Figure |3] displays the binary fraction as a function of 
eccentricity for low inclination objects in the 3:2 and 2:1 
resonances, normalized to the binary fraction of cold clas- 
sical KBOs at 45 AU. These relative fractions are calcu- 
lated by combining the expected formation locations of 
the resonant KBOs from Equation ^ with the trend in 
binary formation rate as a function of formation location 
given by Equation (fT3| and 7 = — 1. The difference be- 
tween true and expected formation location is discussed 
in Section [221 

In general, the larger the semi-major axis of a reso- 
nance, the higher the binary fraction of its cold compo- 
nent should be. To zeroth order, we expect the binary 
fraction in the 2:1 resonance to be roughly comparable 
to the cold classicals. A higher fraction of large eccen- 
tricity objects in resonance would lower this estimate, as 
can be seen from Figure |3l To correct for this bias in a 
given observed sample, one could, for example, calculate 
the estimated fraction using the eccentricity of each ob- 
served object from Figure [5] and calculate the average. 
In contrast, we expect cold 3:2 objects to have a 20-30% 
lower binary fraction than the cold classicals. Again, for 
a particular set of observations, one can calculate an ex- 
pected fraction as above. 

Uncertainties in the population statistics of resonant 
KBOs prevent us from making a robust estimate of the 
number of measurements required to distinguish this ef- 
fect. Nevertheless, we provide the likely scale of the ob- 
servations required by calculating the following example. 
If the binary fraction of cold classicals at 45 AU is ~30%, 
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consiste nt with the observ ed binary fraction among cold 
CKBOs ()Noll et aLll2008bO . then averaged over all eccen- 
tricities, the low inclination 3:2 objects are predicted to 
have binary fraction ~0.22, compared with ~0.27 in the 
2:1 resonance and ^0.29 in the cold classical belt. Iden- 
tifying such subtle differences would require measuring 
the fra ction of eac h population to within approximately 
2-3%. INoll et al.l feOOSb.) measure a binary fraction of 
29^g% in the cold classicals based on binary detection for 
17 of 58 targets. We estimate that a sample containing 
approximately N ~ 250 objects would generate errors of 
y/0.3N /N = 0.035, enough to distinguish our predicted 
binary fractions for cold 3:2 objects from the cold clas- 
sicals. This approximate number of observations would 
be required in each population. The binary fractions in 
the cold 3:2 and cold 2:1 resonant populations could be 
separated with samples of ^-^500 objects each. 

Currently, only 101 3:2 objects are known (as defined 
in ji3.2p . 54 of which have i < 10° and 26 of which have 
i < 5°. In the 2:1 resonance, 15 objects are known, 9 of 
which have i < 1 0° and 6 of whi c h hav e z < 5°. How- 
ever, according to lKavelaars et all (|2009[ ). approximately 
2.51^ 2 ^ 10'' Plutinos are present in the belt with solar 
system absolute magnitude in the g' filter of Hg < 10. 
Given 



D(km) = 2 X 1.5 X l0*^p-"-5l0"-2('"e- 



H) 



(14) 



(e.g.. iPetit et aLl l2008al) with 



7710 



. „ - V " — - -26.47 in g' 

Kavelaars et al T I2009D . H = Hg < 10 corresponds to 
diameters _D > 50 km. We have used albed o p = 0.1, con- 
sisten t with observed KBOs in the visible (jBrucker et al.l 
l2009t) . 

We are interested in objects near the break in the 
KBO size distribution, corresponding to the largest num- 
ber of KBOs that are not likely to have collisionally 
evolved over the age of the solar system. This break 
is measured at app arent magnitude 25 in _R-band 
(jFuentes et al.l [20091 ) . At a distance d from the Sun, an 
object with absolute magnitude H has apparent magni- 
tude H + 2.5 log(d2 AVAU"), with A « d - 1 AU. Using 
d = 42 AU and niQ — —27.12 (appropriate for the Bessel 
R filteiEl) , the break is at D ^ 60km. This differs from 
the 90 km diameter quoted in iFuentes et al.l (|2009t ) due 
to our larger assumed albedo. The actual diameter at 
the break is uncertain to within approximately a factor 
of 2. Since D ~ 50km equals the size at the break to 
within the errors, we conservatively adopt a number of 
Plutinos equal to half of the estimated number of Pluti- 
nos with Hg < 10. Given these assumptions, if a cold 
population of 3:2 KBOs is shown to exist with i < 10° 
and if the fraction of Plutinos with i < 10° is ^^0.5 as in 
the currently observed sample, then the 3:2 population 
likely contains at least ~6 x 10'^ dynamically cold objects 
larger than the break in the size distribution. More than 
an order of magnitude more Plutinos are likely present 
in the belt than would be required to compare the binary 
fractions of cold 3:2 objects and cold classicals. 

The number of objects in the 2:1 resonance is less 
well constrained than th e number in the 3:2 resonance. 
iChiang &: JordanI (|2002l ) estimate that the 2:1 hosts a 
factor of ~3 fewer objects than the 3:2, and Lawler et 

* |http: //mips . as ■ arlzona. edu/~cnaw/suii.html"] August 2010. 



al. (in prep.) estimate a factor of ~4 fewer objects in 
the 2:1 than in the 3:2 based on the well-characterized 
Canada-France Ecliptic Plane Survey. Given a factor of 
4 reduction and assuming, in analogy with the 3:2, that 
half of the 2:1 objects are members of a cold popula- 
tion, a factor of ~3 more 2:1 objects are present than we 
predict are required to differentiate between the binary 
fractions of the cold 3:2 and 2:1 populations. 

LSST is currently projected to detect more than 30,000 
KBOs with diameters larger than ^100 km (I vezic et al.l 
12001 . Using m© = -27.08 (SDSS r' fihci0) and albedo 
p = 0.1, the LSST limiting magnitude of ^ 24.5 in r- 
band corresponds to KBOs with diameters D > 70 km 
at the semi-major axis of the 3:2 resonance {d « 39 AU) 
and D > 100 km at the semi-major axis of the 2:1 reso- 
nance (d w 48 AU). In practice, smaller resonant KBOs 
will be observable because the large eccentricities of these 
objects bring them close r to the Sun. W ith an angu- 
lar resolution of ~0.7" (jlvezic et al.l[2008h . LSST wih 
likely be able to distinguish some wide-separation bina- 
ries without req uirii iR follow-up observ ations . However, 
the observations in iNoll et all ()2008b ') using the Hub- 
ble Space Telescope have a factor of 10 better resolution 
and include binaries t hat will not be resolved by LSST. 
iKavelaars et al.l (|2009l ) estimate that Plutinos are ~20% 
as abundant as CKBOs. If only ~2% or more of the ob- 
jects discovered by LSST are Plutinos and if the fraction 
of wide binaries among these objects is measured, then 
LSST will detect sufficiently many low inclination Pluti- 
nos to search for our predicted difference between the 
binary fractions of cold 3:2 and cold classical objects. 

LSST may also detect the ^-^500 low inclination 2:1 ob- 
jects required to differentiate between the 3:2 and 2:1 
populations. Conservatively choosing a size distribution 
dN/dR (X R"'^ with q — 5 {a.t the large end of the 
q = 2-5 range measured for various subpopulations of 
the belt. IFuentes et al.l 120091: iFraser fc KavelaarsI l200l 
iFraser et al.ll2010[ ). we estimate the number of cold 2:1 
objects with D ^ 100km to be ^200. These objects will 
be observable by LSST at their semi-major axis. A KBO 
in 2:1 resonance with e = 0.2 spends ~16% of its time 
at distances from the Sun of no more than 40 AU. At 
these distances, objects with D ~ 70 km will be visi- 
ble, adding another ~ 800 x 0.16 ~ 100 observable low- 
inclination bodies, for a total of ^300 objects. Given the 
substantial uncertainties encompassed in this calculation 
and our conservative choice of g = 5, it is plausible that 
LSST will detect sufficient low inclination 2:1 bodies to 
differentiate between the cold 3:2 and 2:1 populations. 

We now turn to the search for trends in binary fraction 
with eccentricity within the low inclination population of 
a single resonance. As illustrated in Figure |31 within a 
single resonance, the binary fraction at small eccentric- 
ities is ^--^40% larger than that at the highest observed 
eccentricities. If we consider the average of all objects 
with e < 0.2 compared with those having e > 0.2, the 
difference is reduced to ^-^20%, though this value depends 
on the distribution of eccentricities included in the aver- 
age. 

Extending our example calculation above, we estimate 
that low inclination 3:2 objects have binary fraction 
-0.25 for e < 0.2 and ~0.21 for e > 0.2. These frac- 
tions could be separated with samples of — 750 in 
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each of the e < 0.2 and e > 0.2 subsets of the low- 
inchnation data. Again, these numbers depend on the 
distribution of eccentricities of the objects in the reso- 
nances. Though the observational effort needed would 
be substantial, this experiment should be possible using 
KBOs discovered by LSST. An order of magnitude more 
objects are likely present in the belt than would be re- 
quired to test this effect, and LSST will detect a sufficient 
number of Plutinos to search for our predicted trend if 
^10% or more of its projected KBO discoveries are in 
the 3:2 resonance. 

In our example, the low inclination 2:1 objects are pre- 
dicted to have binary fraction ~0.31 for e < 0.2 and 
^0.26 for e > 0.2. Separation of these fractions re- 
quires '-^lOOO low-inclination 2:1 objects, equally divided 
between low and high eccentricities. We estimate that 
enough objects exist in the belt to perform this exper- 
iment. Though our estimate indicates that LSST will 
discover a factor of ~3 fewer 2:1 KBOs than required, 
we emphasize that this calculation employs a conserva- 
tive size distribution and is highly uncertain because the 
intrinsic population of 100 km objects in the 2:1 reso- 
nance is not well known. 

We recall that because, in the migration scenario, the 
2:1 resonance passed through the cold classical belt, the 
2:1 resonance should provide the clearest case for our dy- 
namical signature. If the 3:2 resonance has a cold compo- 
nent, then it should also provide a good case for testing 
the binary abundance. As mentioned above, higher order 
resonances are less likely to capture unexcited objects. 
However, if they have some cold component, then they 
should show a similar trend in the binary abundance as 
a function of eccentricity as the 3:2 and 2:1 resonances. 

We re-emphasize that these trends should only ex- 
ist among an observed cold (low-inclination) compo- 
nent of the resonant population, if such a component 
exists. We further emphasize that we propose to test 
the binary fraction and not the absolute number of bi- 
naries in the resonant populations. Long-term differ- 
ences in the loss rates from different resonances or as 
a function of orbital parameters within a ny given res- 
onance (e.g., iTiscareno k Malhotral [20091 ) will not af- 
fect our results because these losses should affect bina- 
ries and single KBOs equally. Finally, we note that the 
KBOs need to be tightly in resonance (small libration an- 
gle) rather than tenuously in resonance, since they may 
otherwise have been tran siently captured by scattering 
(jLvkawka fc Mukail 120061 ) . 

5.3. Interpretation for Different Dynamical Sculpting 

Models 

If a cold component is observed in the fraction of wide- 
separation, comparable mass binaries in one or more res- 
onances, it will provide support for migration-induced 
capture, which naturally produces a cold component 
within resonances. Migration-induced capture may be 
consistent with no cold component in the 3:2 resonance, 
but the lack of a cold component in the 2:1 resonance 
would argue strongly against this model. In contrast, i n 
chaotic capture as envisioned bv iLevison et al.l (|2008l ). 
the binary fraction in the resonances should match the 
fraction in the scattered disk, even at low inclinations. 
Given the presence of a distinct binary fraction in the 
cold classical disk, for chaotic capture to be viable, the 



models presented in ILevison et al.l ()2008[ ) need to be ad- 
justed to allow maintenance of a cold classical popula- 
tion that formed in situ. Even if such a modification is 
possible, it would not generate a cold component within 
the resonances, retaining the search for a cold resonant 
component as a useful distinguishing test. 

Measurement of a trend in binary fraction with eccen- 
tricity within the cold component of a resonance would 
provide conclusive evidence of migration-induced cap- 
ture. Unfortunately, lack of such a trend would not con- 
stitute conclusive evidence that migration-induced cap- 
ture did not occur. We expect such a trend to be a 
subtle effect, and uncertainty in the surface-density pro- 
file of the initial disk means that we cannot conclusively 
predict its magnitude. Lack of an observed trend could 
imply a steeper surface density (nearer 7 = —2) or a dif- 
ferent binary formation scenario. Similarly, if a trend in 
binary fraction with eccentricity is observed among low- 
inclination resonant KBOs, but the trend differs from 
our prediction, this could imply either a different surface 
density profile or an alternative binary formation mech- 
anism. 

6. KBO COLORS AND SIZES 

Binary fraction is not the only physical characteristic 
that has been shown to differ between dynamical pop- 
ulations in the Kuiper belt — striking correlations with 
KB O inclination have also been measured for KBO col- 
ors (Tcglcr & Romaiiishinl [20001 
[Peixinho ct al. 20a 



Truiillo & Brown 2002] 
and sizes ([Lcvison & Stern 2001). 
These characteristics are more difficult to interpret than 
the binary fraction because no quantitative theory ex- 
ists for how they should vary as a function of formation 
location. Nevertheless, we now consider the procedure 
outlined in Section [S] in light of these properties. We ad- 
dress the 3:2 resonance, where data is most abundant. A 
detailed reassessment of the available data is warranted 
but beyond the scope of this paper. For our purposes, 
we demonstrate that currently measured size and c olor 
distributions, as c ompiled by ([A lmeida e t al.l I2OO90 for 
3:2 objects and by (jPeixinho et a l. 2008) for classical ob- 
jects, do not conclusively distinguish between our two 
dynamical scenarios. 

First, does the 3:2 resonance contain a cold com- 
ponent? Figures [4| (a) and (b) display KBO diame- 
ters and B-R colors as a function of inclination for ob- 
jects in the 3:2 resonance and for classical objects. We 
convert inclinations to be with respect to the invari- 
able plane. We calculate diameters froi n the absolute 
solar systen i magnitudes Hp quo ted in lAlmeida et al.l 
(|2009f l and iPeixinho et~sl\ (|2008| ) using Equation ^ 
with H = Hr and tuq = -27.12 (Bessel R filtei[3). 
In keeping with the high albed os observed for KBOs in 
the visible ([Brucker et al.l[2009t ). we use al bedo p = 0.1. 
Members of the Haumea collisio n al family (iBrown et al.l 
120071: iRagozzine fc Brownl [20071: iSnodgrass et al.l |2010[) 
are marked. 

Figure m (a) suggests that 3:2 objects with inclinations 
< 5-10° may have systematically smaller sizes than high 
inclination population, consistent with the trend seen 
among classical objects. At first glance, the color trend 
is less promising. However, as noted by lAlmeida et al.l 
(|2009f) . the smallest 3:2 resonant KBOs have systemati- 
cally bluer colors than their larger compatriots (see Fig- 
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Fig. 4. — Observed diameter as a function of inclination with 
respect to the invariable plane (a), inclination vs. B-R color (b), 
and diameter vs. B-R color (c) for classical (gray squares) and 
3:2 resonant (black circles) KBOs. Small circles correspond to 3:2 
objects with sizes smaller than the minimum-size classical KBO. 
Classical objects that are known members of the Haumea coUisional 
family (gray asterisks) have unusually blue colors, known to result 
from water ice featur es. Data for resonant objects is compiled by 
IXlmeida et aP ||2009|) and that for classical KBOs is compiled by 
IPeixinho et al.l I I200SI ). See the text for details of our conversion 
from magnitudes, Hn, to diameters. 

ureH] (c)). The smallest classical object in the sample 



considered here has D ~ 110 km, close t o the observed 
break in the classical size distribution ()Fuentes et al.l 
l2009f) . We suggest that this color difference, unlike the 
color variation among large KBOs, may result from a 
collisional break in the population of 3:2 objects that 
is not yet observed among classical KBOs because col- 
ors have not been measured for small enough objects. 
This possibility merits further i nvestigation (Schlichting 
et al., in prep.). lAlmeida et al.l (|2009) alternatively sug- 
gest that this trend may come from collisions between 
Plutinos and Neptune Trojans. If either idea is vali- 
dated, only the colors of those KBOs too large to suffer 
ongoing collisions reflect conditions in the initial nebula 
and only those should be considered when searching for 
a difference in color between in a low-inclination and a 
high-inclination population. Figure |4] (b) is suggestive, 
though certainly not conclusive, that such a difference 
exists, again divided by a characteristic inclination be- 
tween 5 and 10°. Th is is similar to the color break that 
IPeixinho et"all ()2008[ ) find at ~12° in the classical beU. 
We highlight that the transition in physical properties 
could be at ^^10° or even higher rather than at 5° as 
usually quoted, even in the classical belt. 

Next, do 3:2 objects show a trend in size or color 
with eccentricity? No such trend has been obse rved in 
the 3:2 population (jPoressoundiram et al.l [20081 ). How- 
ever, because no quantitative model exists for how large 
KBO colors vary with formation distance, it is not clear 
whether we would expect to see such a trend in colors. 
For example, the difference in color between cold and hot 
classicals might result not from a continuous variation in 
the protoplanetary nebula but rather from critical tran- 
sitions in the abil ity of KBOs to incorpo rate or retain 
different ices fe.g.. lSchaller fc Brownl l2007). In this case, 
the critical transition separating the region of formation 
for hot and cold CKBOs could have been interior to the 
region swept by the 3:2 resonance, leaving no trend with 
eccentricity in the colors of cold resonant objects. If such 
a trend does exist, we would expect higher eccentricity 
objects to have bluer colors. The number of large, low- 
inclination, 3:2 resonant objects is not large enough to 
determine whether this is the case. 

7. SUMMARY 

Observations demonstrate that unlike asteroids, a sub- 
stantial fraction of Kuiper belt objec ts arc comparable- 
mass, wide separation binaries (jNoll et al. 2008a). These 
bodies typically have sizes of ~100 km. Unlike the Pluto- 
Charon system and other large mass-ratio binaries, they 
likely formed by dynamical capture rather than through 
collisions. Because the Hill radius of a planetesimal in- 
creases with distance from the Sun, the efficiency of dy- 
namical capture is higher for binaries that formed farther 
from the Sun. 

The binary fraction in the cold classical Kuiper belt 
is ~29%, compare d with ~10% for hot classical KBOs 
(|Noll et al.ll2008bD . This discrepancy suggests that ob- 
jects in the hot classical belt formed closer to the Sun 
where binary formation was less efficient. The binary 
fraction of the hot belt may also have been reduced 
by binary disruption during scatteri ng interactions with 
Neptune (jParker fc Kavelaarsll2010f ). In this paper, we 
have developed a method that uses the fraction of wide- 
separation, comparable mass binaries in various popu- 
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lations of the Kuiper belt to distinguish between two 
competing frameworks for the dynamical sculpting of the 
belt. If current projections are correct, we expect LSST 
to discover sufficiently many Kuiper belt objects to per- 
form our proposed tests. 

In both dynamical scenarios, during the late stages of 
planet formation, Neptune was located ~10 AU or more 
closer to the Sun than its current semi-major axis of 30 
AU. A dynamical cataclysm excited the the planet's or- 
bit, which was subsequently damped by dynamical fric- 
tion with planetesimal debris. While Neptune's orbit 
was excited, it interacted with and excited the orbits of 
nearby planetesimal debris. In this context, we consider 
the following two scenarios, named for the methods by 
which resonant KBOs are captured. 

1. "Migration- induced capture." Neptune's orbit was 
damped onto an approximately circular, coplanar 
orbit while it still orbited ^10 AU closer to the 
Sun than it does today. It then migrated smoothly 
outward to its cu rrent location as it sca ttered plan- 
etesimal debris (jFernandez fc Ipiri984[ ). capturing 
KBOs into resonance as it proceeded ([Malhotral 
ll993L[T995l ). Each resonance captured objects from 
a distinct region of the protoplanetary disk, and the 
final eccentricity of a resonant KBO is correlated 
with the location of its capture into resonance. 

2. "Chaotic capture." Neptune scattered quickly onto 
an orbit at most a few AU away from its current 
location. The planet's large eccentricity produced 
a chaotic region of overlapping resonances extend- 
ing out to the 2:1 resonance, allowing objects scat- 
tered by Neptune to fill the region of space cur- 
rently containing the Kuiper belt. Neptune's ec- 
centric ity then da mped, leaving objects in reso- 
nance (jLevison et al. 2008). Due to the chaotic na- 
ture of this process, most or all information about 
the formation locations of KBOs was randomized. 

To distinguish between these models, we calculate the 
fraction of wide-separation, roughly equal-mass binaries 
among KBOs with diameter '^lOO km expected for dif- 
ferent dynamical populations under the hypothesis of 
migration-induced capture. We assume that KBO bi- 
naries formed via transient binary capture, made perma- 
nent by dy namical friction with small bodies (the Li^s 
mechanism. iGoldreich et all 120021 ). We evaluate our ex- 
pressions for binary formation in a disk with surface 
density profile a cx a~^, consistent with observations 
of the outer regions of extrasolar protoplanetary disks 
([Andrews et al.ll2Q09f) . Other binary formation mecha- 
nisms that vary with location in the protoplanetary disk 
and other disk surface density profiles will produce sim- 



ilar, but quantitatively different, signatures. In particu- 
lar, the mechanism produces an opposite trend in the 
binary fraction with formation location compared to the 
L^s mechanism. 

Figure |3] summarizes the pattern of binary fractions 
predicted as a result of migration-induced capture, given 
these choices. A low-inclination component of the 2:1 
resonance should exist with a larger binary fraction than 
its high inclination component. Such a component may 
also exist in the 3:2 resonance. Among bodies with incli- 
nations less than 5-10°, the fraction of binaries in the 2:1 
resonance should be comparable to that in the classical 
belt, while the low-inclination component in the 3:2, if it 
exists, should have a fraction ~20-30% lower than among 
classical KBOs. Within the low inclination component 
of the 3:2 or 2:1 resonance, objects with e < 0.2 should 
have a binary fraction ~20% higher than those with 
e > 0.2. Our calculations are not affected by dynamical 
loss processes that operate with v arying efficiency in dif- 
ferent resonant populations fe.g.. iTiscareno fc Malhotral 
[200i)— these processes affect single and binary KBOs 
equally and we are interested only in the binary fraction. 
Chaotic capture does not produce these signatures. 

Though measurement of the full range of dynamical 
signatures predicted in this paper will require a sub- 
stantial observational effort, searches for low-inclination 
components in the binary fractions of the 3:2 and 2:1 
resonances are within immediate reach. When compared 
with the cold component observed in the classical belt, 
these cold resonant populations need not include the 
same fraction of the overall population, nor must they 
have the same characteristic inclination. Lack of a low 
inclination component in the 2:1 resonance would indi- 
cate that migration-induced capture is unlikely. Con- 
versely, presence of a low inclination component in ei- 
ther resonance would constitute strong evidence for the 
migration-induced capture mechanism. 
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